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1 Introduction 

In a recent series of papers, Nourdin and Peccati [131 114] , and Nourdin, Peccati and Reveillac 
|15j . have shown that one can effectively combine Malliavin calculus on a Gaussian space (see 

■ e.g. [16]) and Stein's method (see e.g. [31 [26]) in order to obtain explicit bounds for the normal 
and non-normal approximation of smooth functionals of Gaussian fields. 

The aim of the present paper is to extend the analysis initiated in [13] to the framework 
of the normal approximation (in the Wasserstein distance) of regular functionals of Poisson 
measures defined on abstract Borel spaces. As in the Gaussian case, the main ingredients of 
^ ■ our analysis are the following: 

CO 

1. A set of Stein differential equations, relating the normal approximation in the Wasser- 
^ \ stein distance to first order differential operators. 

l> 

O I 2. A (Hilbert space- valued) derivative operator D, acting on real- valued square-integrable 

■ random variables. 



operators. 



3. An integration by parts formula, involving the adjoint operator of D. 

4. A "pathwise representation" of D which, in the Poisson case, involves standard difference 



As a by-product of our analysis, we obtain substantial generalizations of the Central Limit 
Theorems (CLTs) for functionals of Poisson measures (for instance, for sequences of single and 
double Wiener-Ito integrals) recently proved in [2H [22j [23] (see also [H [20] for applications 
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of these results to Bayesian non-parametric statistics). In particular, one of the main results 
of the present paper (see Theorem 15.11 below) is a CLT for sequences of multiple Wiener-Ito 
integrals of arbitrary (fixed) order with respect to a general Poisson measure. Our conditions 
are expressed in terms of contraction operators, and can be seen as a Poisson counterpart to 
the CLTs on Wiener space proved by Nualart and Peccati [18] and Nualart and Ortiz-Latorre 
|17j . The reader is referred to Decreusefond and Savy [5] for other applications of Stein- type 
techniques and Malliavin operators to the assessment of Rubinstein distances on configuration 
spaces. 

The remainder of the paper is organized as follows. In Section [21 we discuss some pre- 
liminaries, involving multiplication formulae, Malliavin operators, limit theorems and Stein's 
method. In Section [3l we derive a general inequality concerning the Gaussian approximation 
of regular functionals of Poisson measures. Section H] is devoted to upper bounds for the 
Wasserstein distance, and Section to CLTs for multiple Wiener-Ito integrals of arbitrary 
order. Section [6] deals with sums of a single and a double integral. In Section d we apply our 
results to non-linear functionals of Ornstein-Uhlenbeck Levy processes. 

2 Preliminaries 
2.1 Poisson measures 

Throughout the paper, (Z, Z, fi) indicates a measure space such that Z is a Borel space and \x 
is a cr-finite non-atomic Borel measure. We define the class Z^ as 2^ = {B £ Z : fi (B) < oo}. 
The symbol N = {N (B) : B £ Z^} indicates a compensated Poisson random measure on 
(Z, Z) with control [i. This means that TV is a collection of random variables defined on 
some probability space (Q, J 7 , P), indexed by the elements of 2^, and such that: (i) for every 
B,C G 2^ such that B n C = 0, N (B) and N (C) are independent, (ii) for every B £ 2^, 



where (B) is a Poisson random variable with parameter [J,(B). Note that properties (i)-(ii) 
imply, in particular, that N is an independently scattered (or completely random) measure 
(see e.g. [23J). 

Remark 2.1 According e.g. to [25j Section 1], and due to the assumptions on the space 
(Z,Z,n), it is always possible to define (0, P) in such a way that 



where 5 Z denotes the Dirac mass at z, and N is the compensated canonical mapping given by 



Also, in this case one can take T to be the P-completion of the cr-field generated by the 
mapping (12. 2\ . Note that, under these assumptions, one has that 



P{w : u(B) < oo, \/B s.t. n(B) < oo} = 1 and F{lo : 3z s.t. uj{{z}) > 1} = 0. (2.3) 



N(B) l =<$(B)-n.(B) 




(2.1) 



N(B)(lu)=lu{B)-u,{B), BeZ^. 



(2.2) 
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Remark 2.2 Prom now on, and for the rest of the paper, the hypotheses (|2.2p and (|2.3j) . on 
the structure of Q and N, are implicitly satisfied. In particular, T is the P-completion of the 
a- field generated by the mapping in (|2.2p . 



For every deterministic function h G I? (Z, Z, n) = L 2 (u) we write N (h) = f z h (z) N (dz) 
to indicate the Wiener-Ito integral of h with respect to N (see e.g. [9] or [27] )• We recall 
that, for every h G L 2 (a.), the random variable N (h) has an infinitely divisible law, with 
Levy-Khinchine exponent (see again [27]) given by 

^(M)=logE e iA7?(h) = J e iXh{z) - 1 - iXh (z) fi(dz), Ael. (2.4) 

Recall also the isometric relation: for every g,h G L 2 (//), E[7V (5) iV (/i)] = f z h (z) g (z) \i (dz). 

tix n > 2. We denote by L 2 (/i n ) the space of real valued functions on Z n that are 
square-integrable with respect to /x n , and we write L 2 (fj, n ) to indicate the subspace of L 2 (/j, n ) 
composed of symmetric functions. For every / £ L 2 (/i n ), we denote by / n (/) the multiple 
Wiener-Ito integral of order n, of / with respect to N. Observe that, for every m, n ^ 2, 
/ G L 2 (/i n ) and g G L 2 (/i m ), one has the isometric formula (see e.g. [3T]): 

E [in (/Km (#)] = ra!(/,5) L 2 (At n)l (n=m) . (2.5) 

Now fix n ^ 2 and / G L 2 (fj, n ) (not necessarily symmetric) and denote by / the canonical 
symmetrization of /: for future use we stress that, as a consequence of Jensen inequality, 

Wfh^) < II/IU^)- (2-6) 

The Hilbert space of random variables of the type I n (f), where n ^ 1 and / G L 2 (/i n ) is called 
the nth Wiener chaos associated with N. We also use the following conventional notation: 
h(f) = N(f), f G L 2 (fi); I n {f) = /„,(/), / G L 2 (fi n ), n > 2 (this convention extends 
the definition of I n {f) to non-symmetric functions /); Io(c) = c, c G M. The following 
proposition, whose content is known as the chaotic representation property of N, is one of the 
crucial results used in this paper. See e.g. [TH] or [32] for a proof. 

Proposition 2.3 (Chaotic decomposition) Every ran dom variable F G L 2 (T,F) = L 2 (P) 
admits a (unique) chaotic decomposition of the type 



F = E(F) + J2 I n(fn), (2-7) 
where the series converges in L 2 and, for each n ^ 1, the kernel f n is an element of L 2 (/i n ). 



2.2 Contractions, stars and products 

We now recall a useful version of the multiplication formula for multiple Poisson integrals. To 
this end, we define, for q,p ^ 1, / G L 2 (pi p ), g G L 2 (u 9 ), r = 0, ...,qAp and Z = l,...,r, 
the (contraction) kernel on Z p+q ~ r ~\ which reduces the number of variables in the product 
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fg from p + q to p + q — r — I as follows: r variables are identified and, among these, I are 
integrated out. This contraction kernel is formally defined as follows: 

f* r 3(7li • • • Hr-hh,-- ■ ,tp-r,Sl, . . . ,Sg- r ) (2.8) 

= / f{z 1 ,...,zi,^ 1 ,...,^ r -i,t 1 ,...,tp- r )g(z 1 ,...,zi,^ 1 ,...,^ r _i,s 1 ,...,s q -r)^(dz 1 ...dzi), 
Jz' 

and, for I = 0, 

/*rff(7l) • • • ,7n*l, • • • ,tp-r,Sl, • • • , Sg-r) = /(7l, • • • i7n*li • • • ,tp-r)g{ll, ■ ■ ■ ,7r,Sl, ■ • • 

(2.9) 

so that /*g (/(ti, . . . ,i p ,si, . . . ,s 9 ) = /(ti, . . . ,t p )g(s 1} . . .,s q ). For example, if p = q = 2, 

/ *° 5 (7, M) = / (7, *) 5 (7, s), f*\g(t,s) = [ f (z, t) g (z, s) p (dz) (2.10) 

f4g(l)= [ f(z, 1 )g(z, 1 )^(dz) (2.11) 
Jz 

f*29= / f(zi,Z2)g(z 1 ,Z2)f-i(dzi)n(dz 2 ). (2.12) 



z Jz 



The quite suggestive 'star-type' notation is standard, and has been first used by Kabanov in 
[8] (but see also Surgailis [31]). Plainly, for some choice of f,g,r,l the contraction f * l r g may 
not exist, in the sense that its definition involves integrals that are not well-defined. On the 
positive side, the contractions of the following three types are well-defined (although possibly 
infinite) for every q ^ 2 and every kernel / G L 2 (/j, q ): 

(a) / *° f(zi, Z2q- r ), where r = 0, q, as obtained from (|2.9p . by setting g = /; 

(b) f* q f(zi, ...,z q -i) = J z i f 2 (zi, ...,Zq-i,-)dn l , for every I = 1, ...,q; 

(c) /*£ /, for r = l,....,q- 1. 

In particular, a contraction of the type / *_ /, where / = 1, ...,q — 1 may equal +oo at some 
point (z\, z q -i). The following (elementary) statement ensures that any kernel of the type 
/ g is square-integrable. 

Lemma 2.4 Let p,q ^ 1, and /ei / G L 2 (n q ) and g G L 2 (/i p ). Fix r = 0,...,q A p. Then, 
f* r r g£L 2 (^+i- 2r ). 



Proof. Just use equation (|2.8p in the case I = r, and deduce the conclusion by a standard 
use of the Cauchy-Schwarz inequality. ^ 

We also record the following identity (which is easily verified by a standard Fubini argu- 
ment), valid for every q ^ 1, every p = 1, q: 

(f4f) 2 dvP>-*= f (/*«-* /)V, (2.i3) 

Z 2 i-p JZP 

for every / G L 2 (fj, q ). The forthcoming product formula for two Poisson multiple integrals is 
proved e.g. in [8] and |31j . 
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Proposition 2.5 (Product formula) Let f G L 2 (fi p ) and g G L 2 (/i 9 ), p, q ^ 1, and sup- 
pose moreover that f -k r g G L 2 (/jL p+q ~ r ~ l ) for every r = 0, A g and / = 1, r s?«c/i i/ia£ 
Z 7^ r. T/ien, 

= (*) e (;)vh(/Ts), (2.i4) 

where the tilde "~ " stands for symmetrization, that is, 

f * l r g (xi, ...,X q+p _ r _i) = ^ q l p ~ fj\^2f* l r 9 {Xa(l),-,Xq( q+ p-r-l)) , 

where a runs over all (q + p — r — 1)1 permutations of the set {1, q + p — r — I}. 

Remark 2.6 (Multiple integrals and Levy processes). In the multiple Wiener-Ito integrals 
introduced in this section, the integrators are compensated Poisson measures of the type 
N(dz), defined on some abstract Borel space. It is well-known that one can also build similar 
objects in the framework of Levy processes indexed by the real line. Suppose indeed that we 
are given a cadlag Levy process X = {X t , t > 0} (that is, X has stationary and independent 
increments, X is continuous in probability and Xq =0), defined on a complete probability 
space (0, P), and with Levy triplet (7, a 2 , v) where 7 G M, a > and v is a Levy measure 
on R (see e.g. [27]). The process X admits a Levy-Ito representation 

X t = 1 t + aW t + [ xdN(s,x) +lim // xdN(s,x), (2.15) 

J J (0,t]x{\x\>l} £ i° J J (0,t]x{e<\x\<l} 

where: (i) {Wt, t > 0} is a standard Brownian motion, (ii) 

N(B) = #{t : (t, AX t ) G B}, Be B((0, 00) x M ), 

is the jump measure associated with X (where Mo = R — {0}, AX t = X t — X t - and j^A 
denotes the cardinal of a set A), and (iii) 

dN(t, x) = dN(t, x) - dt dv(x) 

is the compensated jump measure. The convergence in (|2.15p is a.s. uniform (in t) on every 
bounded interval. Following Ito [7], the process X can be extended to an independently 
scattered random measure M on (M + x R, B(R + x R)) as follows. First, consider the measure 
dfi*(t, x) = o~ 2 dtddo(x) + x 2 dtdv(x), where 5q is the Dirac measure at point 0, and dt is the 
Lebesgue measure on R. This means that, for E G £>(R + x R), 

fi*(E)=a 2 [ dt+ [f x 2 dtdv(x), 
Je(o) JJe> 

where E(0) = {t G R+ : (i,0) G E} and E' = E - {(t,0) G E}; this measure is continuous 
(see Ito [7], p. 256). Now, for E G £(R+ x R) with (j*(E) < 00, define 

M(E) = a f dW t + lim ff x dN(t, x), 

Je(0) n J J{(t,x)eE:l/n<\x\<n} 
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(with convergence in L 2 (Q)), that is, M is a centered independently scattered random measure 
such that 

E[M(E 1 )M(E 2 )] = n E 2 ), 

for E 1 ,E 2 G B(R + x R) with n*(Ei) < oo and /i*(E 2 ) < oo. Now write L 2 n = L 2 {(R + x 
M) n ,,6(IR + x R) n , / u*® n ). For every / G L 2 , one can now define a multiple stochastic integral 
In(f), with integrator M, by using the same procedure as in the Poisson case (the first 
construction of this type is due to Ito - see [7]). Since M is defined on a product space, in this 
scenario we can separate the time and the jump sizes. If the Levy process has no Brownian 
components, then M has the representation M(dz) = M(ds,dx) = xN(ds,dx), that is, M 
is obtained by integrating a factor x with respect to the underlying compensated Poisson 
measure of intensity dtv(dx). Therefore, in order to define multiple integrals of order n with 
respect to M, one should consider kernels that are square-integrable with respects the measure 
(x 2 dt v(dx))® n . The product formula of multiple integrals in this different context would be 
analogous to the one given in Proposition 12.51 but in the definitions of the contraction kernels 
one has to to take into account the supplementary factor x (see e.g. [10] for more details on 
this point). 



2.3 Four Malliavin-type operators 

In this section, we introduce four operators of the Malliavin-type, that are involved in the 
estimates of the subsequent sections. Each of these operators is defined in terms of the chaotic 
expansions of the elements in its domain, thus implicitly exploiting the fact that the chaotic 
representation property (|2.7|) induces an isomorphism between L 2 (P) and the symmetric Fock 
space canonically associated with L 2 (/i). The reader is referred to [19] or [28] for more 
details on the construction of these operators, as well as for further relations with chaotic 
expansions and Fock spaces. One crucial fact in our analysis is that the derivative operator D 
(to be formally introduced below) admits a neat characterization in terms of a usual difference 
operator (see the forthcoming Lemma l2.7p . 

For later reference, we recall that the space L 2 (P;L 2 (^l)) ~ L 2 (Q x Z, JF® Z, P(g> //) is the 
space of the measurable random functions m:1!xZ->1 such that 



E 



u 2 \i{dz] 



z 



< oo. 



In what follows, given / G L 2 (fi q ) (q ^ 2) and z G Z, we write f(z, •) to indicate the function 
on Z q ~ l given by (zi, z,_i) -> f(z, zi, 

i) The derivative operator D. The derivative operator, denoted by D, transforms random 
variables into random functions. Formally, the domain of D, written domD, is the set of those 
random variables F G L 2 (P) admitting a chaotic decomposition (|2.7j) such that 

E nn! n/"HW») <o °- ( 2 - 16 ) 

If F verifies (12.161) (that is, if F G domD), then the random function z — ► D Z F is given by 
D z F = J2nI n -i{f{z,-)), zeZ. (2.17) 
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For instance, if F = ii(/), then D Z F is the non-random function z — ► f(z). If F = h(f), 
then D Z F is the random function 

z-2ii(/(z,.)). (2.18) 

By exploiting the isometric properties of multiple integrals, and thanks to (|2.16p . one sees 
immediately that DF G L 2 (P;L 2 (^)), for every F G domD. Now recall that, in this paper, 
the underlying probability space (O, J 7 , P) is such that O is the collection of discrete measures 
given by (|2.1|) . Fix z£ Z; given a random variable F : — > R, we define F z to be the random 
variable obtained by adding the Dirac mass 5 Z to the argument of F. Formally, for every 
uj G 0, we set 

F z (lu) = F{u + 5 z ). (2.19) 

A version of the following result, whose content is one of the staples of our analysis, is proved 
e.g. in [HI Theorem 6.2] (see also [2H Proposition 5.1]). One should also note that the proof 
given in |19| applies to the framework of a Radon control measure /i: however, the arguments 
extend immediately to the case of a Borel control measure considered in the present paper, 
and we do not reproduce them here (see, however, the remarks below). It relates F z to D Z F 
and provides a representation of D as a difference operator. 

Lemma 2.7 For every F G domD, 

D Z F = F z -F, a.e.-n(dz). (2.20) 



Remark 2.8 1. In the Ito-type framework detailed in Remark [2.61 we could alternatively 
use the definition of the Malliavin derivative introduced in [28], that is, 

Fn x \ — F „ 

D(f >x )F = — 1 , a.e.—x dtv(dx). 

See also |llj . 

2. When applied to the case F = h(f), formula (|2.20|) is just a consequence of the straight- 
forward relation 

F z -F= [ f(x)(N(dx) + S z (dx)) - [ f(x)N(dx) = [ f(x)5 z (dx) = f(z). 
Jz Jz Jz 

3. To have an intuition of the proof of (|2.20|) in the general case, consider for instance a 
random variable of the type F = N(A) x N(B), where the sets A, B G are disjoint. 
Then, one has that F = 12(f), where f(x,y) = 2~ 1 {1a(x)1b(u) + 1a(2/)1_b(^)}- Using 
([2TT8D . we have that 

D Z F = N(B)l A (z) + N(A)l B (z), z G Z, (2.21) 

and 

F Z -F = {N + S Z }(A) x{N + 5 Z }(B) - N(A)N(B), z G Z. 
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Now fix z G Z . There are three possible cases: (i) z £ A and z ^ B, (ii) z £ A, and (iii) 
z £ B. If 2 is as in (i), then F z = F . If z is as in (ii) (resp. as in (hi)), then 

F z -F = {N + 5 Z }{A) x N(B) - N(A)N(B) = N(B) 

(resp. 

F z -F = N(A) x{N + 5 Z }(B) - N(A)N(B) = N(A) ). 

As a consequence, one deduces that F z — F = N(B)1a(z) + N(A)1b(z), so that relation 
(123U1) is obtained from (123TT) . 

4. Observe also that Lemma l2~Tl yields that, if F, G G domD are such that FG G domD, 
then 

D(FG) = FDG + GDF + Z?GL»F, (2.22) 
(see [IH Lemma 6.1] for a detailed proof of this fact). 

ii) The Skorohod integral 5. Observe that, due to the chaotic representation property of TV, 
every random function u G L 2 (P, L 2 (fi)) admits a (unique) representation of the type 



u z = 

n>0 



J^n(/nCv)), ^G^, (2.23) 

where, for every z, the kernel f n (z, ■) is an element of L 2 s ([i n ). The domain of the Skorohod 
integral operator, denoted by domJ, is defined as the collections of those u G L 2 (P, L 2 (/i)) 
such that the chaotic expansion (12, 23ft verifies the condition 

X> + l)!||/„||i 20 , n+1) <oo. (2.24) 
If u G dom<5, then the random variable 6(u) is defined as 

5(u) = J2 I n+l(fn), (2.25) 

where f n stands for the canonical symmetrization of / (as a function in n + 1 variables). 
For instance, if u{z) = f(z) is a deterministic element of £ 2 (/i), then 5(u) = ii(/). If 
u(z) = Ii(f(z,-)), with / G L 2 (fi 2 ), then S(u) = hif) (we stress that we have assumed / to 
be symmetric). The following classic result, proved e.g. in [19], provides a characterization of 
S as the adjoint of the derivative D. 

Lemma 2.9 (Integration by parts formula) For every G G domD and every u G dom<5, 
one has that 

E[G5(u)]=E[(DG,u) L 2 M ], (2.26) 

where 

(DG, u)l 2 {ji) = J D * G x u(z)fj,(dz). 



S 



iii) The Ornstein-Uhlenbeck generator L. The domain of the Ornstein-Uhlenbeck generator 
(see [TBI Chapter 1]), written domL, is given by those F G L 2 (P) such that their chaotic 
expansion (|2.7|) verifies 

^2 n2n] \\fn\\ 2 L 2 {lMn) < OO. 

If F G domL, then the random variable Li* 1 is given by 

LF = -J2 nI n(fn)- (2.27) 

Note that E(LL) = 0, by definition. The following result is a direct consequence of the 
definitions of L>, 5 and L. 

Lemma 2.10 For every F G domL, one has that F G domL) and DF G domJ. Moreover, 

5DF = -LF. (2.28) 

Proof. The first part of the statement is easily proved by applying the definitions of domD 
and dom<5 given above. In view of Proposition 12.31 it is now enough to prove (|2.28|) for 
a random variable of the type F = I q (f), q ^ 1. In this case, one has immediately that 
D Z F = ql q -i(f(z, •)), so that 5DF = ql q (f) = -LF. □ 

iv) The inverse of L . The domain of L , denoted by Lq(P), is the space of centered random 
variables in L 2 (P). If F G Lq(P) and F = Y^n^i In(fn) (and thus is centered) then 

L- l F = - V -/„(/„). (2.29) 
z — ' n 



2.4 Normal approximation in the Wasserstein distance, via Stein's method 

We shall now give a short account of Stein's method, as applied to normal approximations in 
the Wasserstein distance. We denote by Lip(l) the class of real-valued Lipschitz functions, 
from M to M, with Lipschitz constant less or equal to one, that is, functions h that are absolutely 
continuous and satisfy the relation ||fo'||oo ^ 1- Given two real- valued random variables U and 
Y, the Wasserstein distance between the laws of U and Y, written dw(U, Y) is defined as 

d w {U,Y)= sup \E[h(U)]-K[h(Y)]\. (2.30) 
/ieLip(i) 

We recall that the topology induced by dw on the class of probability measures over M is 
strictly stronger than the topology of weak convergence (see e.g. [6]). We will denote by 
,/K(0, 1) the law of a centered standard Gaussian random variable. 

Now let X ~ ^(0, 1). Consider a real- valued function h : M — > R such that the expectation 
E[/i(X)] is well-defined. The Stein equation associated with h and X is classically given by 

h(x)-E[h(X)} = f'(x) - xf{x), xGE. (2.31) 

A solution to (I2.3ip is a function / depending on h which is Lebesgue a.e.-differentiable, and 
such that there exists a version of /' verifying (|2.31|) for every x G P. The following result 
is basically due to Stein [30] . The proof of point (i) (whose content is usually referred as 
Stein's lemma) involves a standard use of the Fubini theorem (see e.g. [3J Lemma 2.1]). Point 
(ii) is proved e.g. in [2j Lemma 4.3]. 
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Lemma 2.11 (i) Let W be a random variable. Then, W = X ~ ^(0, 1) if, and only if, 



E[f(W) - Wf(W)} = 0, 



(2.32) 



for every continuous and piecewise continuously differentiable function f verifying the 
relation E\f'(X)\ < oo. 

(ii) // h is absolutely continuous with bounded derivative, then \2. 31}) has a solution fh which 
is twice differentiable and such that ||/^||oo ^ ll^'lloo and ||/^||oo ^ 2 ] | /i' ] | qq . 

Let J^w denote the class of twice differentiable functions, whose first derivative is bounded 
by 1 and whose second derivative is bounded by 2. If h € Lip(l) and thus ||/i'||oo ^ 1, then 
the solution fh, appearing in Lemma f2. 11( h). satisfies ||/^||oo ^ 1 and ||/^||oo ^ 2, and hence 



Now consider a Gaussian random variable X ~ o/K(0, 1), and let Y be a generic random 
variable such that EY 2 < oo. By integrating both sides of (|2.31|) with respect to the law of Y 
and by using (I2.30P , one sees immediately that Point (ii) of Lemma 12.111 implies that 



Note that the square-integrability of Y implies that the quantity E[y_/*(Y)] is well defined 
(recall that / is Lipschitz). In the subsequent sections, we will show that one can effectively 
bound the quantity appearing on the RHS of (|2.33p by means of the operators introduced in 
Section 12.31 

3 A general inequality involving Poisson functionals 

The following estimate, involving the normal approximation of smooth functionals of N, will 
be crucial for the rest of the paper. We use the notation introduced in the previous section. 

Theorem 3.1 Let F £ domD be such that E(F) = 0. Let X ~ jV(Q, 1). Then, 



fh G &w. 



d w {Y,X)^ sup \E(f'(Y)-Yf(Y))\. 



(2.33) 



fe&w 




where we used the standard notation 





Jz Jz 



(3-3) 



(3.4) 
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Proof. By virtue of the Stein-type bound (|2.33j) . it is sufficient to prove that, for every 
function / such that H/'Hoo ^ 1 and ||/"||oo ^ 2 (that is, for every / S ^w), the quantity 
\E[f(F)-Ff(F)]\ is smaller than the RHS of ([32]). To see this, fix / € J? w and observe that, 
by using (|2.19p . for every to and every z one has that D z f(F){w) = f(F) z (ui) — f(F)(to) = 
f(F z )(u) — f(F)(uj). Now use twice Lemma [2~7l combined with a standard Taylor expansion, 
and write 

D z f(F) = f(F) z -f{F) = f{F z )-f{F) (3.5) 
= f'(F)(F z -F) + R(F Z — F) = f(F){D z F) + R(D Z F), 

where (due to the fact that ||/"||oo ^ 2) the mapping y — > R(y) is such that \R(y)\ ^ y 2 for 
every y £ R. We can therefore apply (in order) Lemma 12.101 and Lemma 12.91 (in the case 
u = DL- l F and G = f(Fj) and infer that 

E[Ff(F)} = E[LL~ 1 Ff(F)} = E[-S(DL- 1 F)f(F)} = E[(Df(F),-DL- 1 F) L 2 M }. 

According to (|3.5p . one has that 

E[(Df(F), -DL^F)^} = E[f(F)(DF, -DL^F)^] + E[(R(DF), -DL^F)^}. 

It follows that 

\E[f(F) - Ff(F)]\ ^ |E[/'(F)(1 - (DF, - D L' 1 F) L2 + \E[(R(DF), -DL^F)^]]. 
By the fact that ||/'||oo ^ 1 and by Cauchy-Schwarz, 

\E[f'(F)(l- (DF,-DL- l F) L2(ll) )]\ < E[|l - [DF,-DL~ l F) L 2^\] 

< [{l-(DF,-DL-iF) L 2 M y]. 
On the other hand, one sees immediately that 

\E[(R(DF),-DL- 1 F) L 2 M }\ < j E [\R(D Z F)D Z L~ 1 F\] fj,(dz) 

Jz 

< / E^D.F^D^FW^dz), 
Jz 

Relations (|3.3p and (|3.4|) are immediate consequences of the definition of L~ l given in (|2.29|) . 
The proof is complete. □ 

Remark 3.2 Note that, in general, the bounds in formulae (I3.1|) - (|3.2p can be infinite. In the 
forthcoming Sections [4HZl we will exhibit several examples of random variables in donxD such 
that the bounds in the statement of Theorem 13.11 are finite. 

Remark 3.3 Let G be an isonormal Gaussian process (see [16]) over some separable Hilbert 
space fj, and suppose that F € L 2 (a(G)) is centered and differentiable in the sense of Malli- 
avin. Then, the Malliavin derivative of F, noted DF, is a i5-valued random element, and in 
[T3] it is proved that one has the following upper bound on the Wasserstein distance between 
the law of F and the law of X ~ o/K(0, 1): 

d w {F,X) <E|1- (DF, —DL^F)^], 

where L~ l is the inverse of the Ornstein-Uhlenbeck generator associated with G. 
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Now fix h G L 2 (fi). It is clear that the random variable I\{h) = N(h) is an element 
of domD. One has that D z N(h) = h(z); moreover -L~ l N{h) = N(h). Thanks to these 
relations, one deduces immediately from Theorem 13.11 the following refinement of Part A of 
Theorem 3 in 1231. 



Corollary 3.4 Let h £ L (/x) and let X ~ ^V(0, 1). Then, the following bound holds: 

d w (N(h),X) < |l- ||/i||i 2w |+y \h(z)\ 3 fi(dz). (3.6) 
As a consequence, if fJ-(Z) = oo and if a sequence {hk} C L 2 (n) n L 3 (p) verifies, as k — > oo, 

IIMl 2 ^) 1 anc? IIMU 3 (aO ^ °> ( 3J ) 

one has the CLT 

N(h k )^X, (3.8) 

and the inequality \3. 6]) provides an explicit upper bound in the Wasserstein distance. 

In the following section, we will use Theorem 13.11 in order to prove general bounds for 
the normal approximation of multiple integrals of arbitrary order. As a preparation, we now 
present two simple applications of Corollary 13.41 

Example 3.5 Consider a centered Poisson measure JV on Z = M + , with control measure p» 
equal to the Lebesgue measure. Then, the random variable fc -1 / 2 iV([0, k)) = N(hk), where 
hk = & -1 ^ 2 l[o.fc]> is an element of the first Wiener chaos associated with N. Plainly, since the 
random variables N((i — (i = l,...,k) are i.i.d. centered Poisson with unitary variance, 

a standard application of the Central Limit Theorem yields that, as k — ► oo, N(hk) ^+ X ~ 
^V(0,1). Moreover, N(hk) G domD for every k, and DN{hk) = hk- Since 



IIMia&O = 1 and j \h k \ 3 fi(dz) = 



one deduces from Corollary 13.41 that 



d w (N(h k ),X) «C 



jfc!/2' 

which is consistent with the usual Berry-Esseen estimates. 



□ 



Example 3.6 Fix A > 0. We consider the Ornstein-Uhlenbeck Levy process given by 

Y t x = V2\[ I uexp(-\(t-x))N(du,dx) , t^O, (3.9) 

where is a centered Poisson measure over RxK, with control measure given by fj,(du, dx) = 
v (du) dx, where u (•) is positive, non-atomic and normalized in such a way that f R u 2 dv = 1. 
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We assume also that J R |u| 3 dv < oo. Note that Y t x is a stationary moving average Levy 
process. According to [231 Theorem 4], on has that, as T — > oo, 



Y t x dt ^ X ~ JT (0, 1) . 



We shall show that Corollary 13.41 implies the existence of a finite constant q\ > 0, depending 
uniquely on A, such that, for every T > 



d w (A x ,X) < 



gA 
T l/2- 



(3.10) 



To see this, first use a Pubini theorem in order to represent Aj, as an integral with respect to 
N, that is, write 



A? 



T 



U 



\2TI\) 



exp (—A (t — x)) dt 



xVO 



N(du,dx) := N(h$). 



Clearly, Aj, S domD and DA^, = hj,. By using the computations contained in [23j Proof of 
Theorem 4], one sees immediately that there exists a constant f3\, depending uniquely on A 
and such that 



/ 



\hT(u,x)\^v(du)dx < 



/9a 

T l/2- 



Since one has also that |||^tIIl 2 (/^) — ^1 = the estimate (I3.10p is immediately deduced 

from Corollary 13.41 □ 



4 Bounds on the Wasserstein distance for multiple integrals of 
arbitrary order 

In this section we establish general upper bounds for multiple Wiener-Ito integrals of arbitrary 
order q ^ 2, with respect to the compensated Poisson measure N. Our techniques hinge on 
the forthcoming Theorem 14.21 which uses the product formula (I2.14p and the inequality (13.21) . 

4.1 The operators G q p and G q p 

Fix q ^ 2 and let / G L 2 s (ii q ). The operator Gp transforms the function /, of q variables, into 
a function Gpf of p variables, where p can be as large as 2q. When p = 0, we set 

G Q f = ql\\f\\l H ^ 

and, for every p = 1, ...,2g, we define the function (z\, ...,z p ) — ► Gpf(z\, ■■■,z p ), from Z p into 
M, as follows: 

Glf( Zl ,...,z p ) = iziZ 1 {^i= p} r\( q ) ( r ^f^f(z 1 ,... : z p ), (4.1) 

r=0 1=0 V / V / 

where the 'star' contractions have been defined in formulae (|2.8p and (|2.9|) . and the tilde "~" 
indicates a symmetrization. The notation (|4.1|) is mainly introduced in order to give a more 
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compact representation of the RHS of (12. 14ft when g = f. Indeed, suppose that / G L 2 (/i 9 ) 
(q ^ 2), and that f * l r f G L 2 (/i 2g ~ r ~') for every r = 0, g and Z = 1, r such that Z / r; 
then, by using (|2.14|) and (14. ip . one deduces that 

W) 2 = E^p(G|/). (4-2) 

p=0 

where I (G 9 Q f) = G q f = g!||/||^ ( 9) . Note that the advantage of (g2J) (over d2~T4l ) is that 
the square I q (f) 2 is now represented as an orthogonal sum of multiple integrals. As before, 
given / G L 2 (n q ) (q ^ 2) and z G Z, we write /(z, •) to indicate the function on given 
by (zi, — > /(z, Zi, z q -i). To simplify the presentation of the forthcoming results, 
we now introduce a further assumption. 

Assumption A. For the rest of the paper, whenever considering a kernel / G L 2 (/j, q ), we will 
implicitly assume that every contraction of the type 

(Zi, ...,Z 2q - r -l) > \ f\ *r \ f\(zi,-,Z2q-r-l) 

is well-defined and finite for every r = l,...,q, every I = l,...,r and every (zi, Z2 q - r -i) G 

Assumption A ensures, in particular, that the following relations are true for every r = 
0, ....,<? - 1, every Z = 0, ...,r and every (zi, z 2 ( 9 -i)-r-0 G Z 2 ( <? ~ 1 )~ r ~': 



[/(z,-)4/(^-)M^) = {/*^ 1 i/} ) and / [/(*,•) 4 /(*r)M<te) = U*r+\/ 



(note that the symmetrization on the LHS of the second equality does not involve the variable 
*)• 

The notation Gp /(z, •) (p = 0, 2(q — 1)) stands for the action of the operator Gp 
(defined according to (14.11) ) on f(z, •), that is, 



G*-V(z, ■)(*!, (4-3) 
S £ 1 {2( 9 -l)-»-l=p} r! r J ^jfiz,-)^ f(zr){zi,-..,Zp). 



r=0 1=0 

For instance, if / G L 2 ^ 3 ), then 



G 2 f(z, -)(a) = 4 x /(z, •) ^ /(z, -)(a) = 4 / /(z, a, u) 2 M (d«). 

Jz 

Note also that, for a fixed z G Z, the quantity Gq _1 /(z, •) is given by the following constant 



Gl 1 f(z,-) = (q-l)l / f\z,-W q -\ (4.4) 
Jz*- 1 

Finally, for every q ^ 2 and every / G L 2 (u. q ), we set 

35/= / GT 1 /(^-H^) = 9- 1 ^/ = (g-l)!||/||i 2(M9) , 
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and, for p = 1, ...,2(q — 1), we define the function (z\, z p ) — ► Gpf(z±, z p ), from Z p into 
R, as follows: 



G q p f(-)= / G p r ^ 1 f(z, -)fi(dz), (4.5) 
Jz 

or, more explicitly, 

Gpf{zi, z p ) 

= /EEVi)- r -M r! r (M/(z,-)4/(^-)(^i,-,%)M^) (4-6) 



r=0 (=0 

g *A(g-l) , 1N 2 



where in (14. 7p we have used the change of variables t = r + 1 and s = I + 1 , as well as the fact 
that p ^ 1. We stress that the symmetrization in (|4.6|) does not involve the variable z. 

4.2 Bounds on the Wasserstein distance 

When q ^ 2 and /x(Z) = oo, we shall focus on kernels / £ L 2 {^jL q ) verifying the following 
technical condition: for every p = 1, 2(g — 1), 



{GT 1 /^,-)} 2 ^ 

Zp 



fi(dz) < oo. (4.8) 



As will become clear from the subsequent discussion, the requirement (|4.8|) is used to justify 
a Fubini argument. 

Remark 4.1 1. When q = 2, one deduces from (|4.3|) that Gj/(z,-)(x) = f(z,x) 2 and 
G\f(z, -)(x, y) = f(z, x)f(z, y). It follows that, in this case, condition (14. 8ft is verified if, 
and only if, the following relation holds: 



z V Jz 



f(z,a) A fi(da) fi(dz) < oo. (4.9) 
Indeed, since / is square-integrable, the additional relation 



f(z,a) 2 f(z,b) 2 fi(da)n(db) fi(dz) < oo (4.10) 

z V Jz 2 
is always satisfied, since 



f(z,a) 2 f(z,b) 2 ^da)fi(db)fi(dz) = ||/||£ 2(/x2) . 
z V 
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2. From relation (|4.3j) . one deduces immediately that (|4.8p is implied by the following 
(stronger) condition: for every p = 1, 2(g— 1) and every (r, Z) such that 2(g— 1) — r — I = 
P 



{f(z,-)4f(z,-)}2 d/J> 



fi(dz) < oo. 



(4.11) 



3. When fjt(Z) < oo and 



/ / {Gfifiz, -)} 2 di? 
jz Uzp 



fx(dz) < oo, 



condition f!4.8|) is automatically satisfied (to see this, just apply the Cauchy-Schwarz 
inequality) . 



4. Arguing as in the previous point, a sufficient condition for (|4.8p to be satisfied, is that 
the support of the symmetric function / is contained in a set of the type A x • • • x A 
where A is such that fJ>(A) < oo . 

Theorem 4.2 (Wasserstein bounds on a fixed chaos) Fix q ^ 2 and let X ~ ^(0, 1). 
Let f G L 2 [ (fi q ) be such that: 



(i) whenever n(Z) = oo, condition (4-8) is satisfied for every p = 1, ...,2(g — 1); 



(ii) for d[i-almost every z £ Z, every r = l,...,q — 1 and every I = 0, ...,r — 1, the kernel 
f(z,-)* l r f(z,-) is an element of L 2 ,(^ 2 ( 9 ~ 1 )~ r ~ i ). 

Denote by I q (f) the multiple Wiener-Ito integral, of order q, of f with respect to N . Then, 
the following bound holds: 



d w (I q (f),X) ^ 



(i-tfii/ii? 



2(9-1) 2 

+ q 2 E W { d pf} d ^ P 

p=l ^ ZP 



2(9-1) f 

EW G q p ~ 1 f{z,-) 2 dpP\ii{dz), 



(4.12) 



(4.13) 



where the notations Gpf and Gp 1 f(z,-) are defined, respectively, in |^.5[ )- |7^1 ) and |^.3|) . 
Moreover, the bound appearing on the RHS of can be assessed by means of the 

following estimate: 



\ 



2(9-1) 



(wiii/ii? 



;=1 ^ p 1 J 



^ II - q\ 



9 *A(g-l) 

+ «E E 1 {2^2,-i}(2g-t- S )! 1/2 (t-l)! 1/2 x 

i=l s=l 

g - IN A - 1 



(4.14) 



(4.15) 



t - 1 J \s - 1 



1/2 



*f /IU 2 ou 2 «-*-») ; 
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Also, if one has that 

/*<f 6 /GL 2 (A V6 = l,...,9-1, (4.16) 

then, 



\ 



2 (9- 1 ) ,■ r r >. 2 



E p! ^ { ^ Gl~ l f{z, -) 2 ^} < (4-17) 

q 6-1 

E E 1 {i^+K2 9 -i}(« + & ) jl/2 (9 - « - !) !l/2 >< (4-18) 



6=1 a=0 



x C-T-.)(^-r) I/> »^""<^>- 



Remark 4.3 1. There are contraction norms in (I4.18P that do not appear in the previous 
formula (|4.15j) . and vice versa. For example, in (14. 18ft one has ||/*° f\\L 2 U 2 i- b )) where 
b = l,...,q (this corresponds to the case b £ {l,...,q} and a = 0). By using formula 
(|2.13l) . these norms can be computed as follows: 

Wf4f\\vQfr->) = \\f*r b f\\W)> b = l,...,q-l; (4.19) 
\\f*° q f\\v>M = J/zJ 4 ^"- ( 42 °) 

We stress that f 2 = f*° q f, and therefore /*° f G L 2 (^) if, and only if, / G L 4 (^). 

2. One should compare Theorem 14.21 with Proposition 3.2 in [13], which provides upper 
bounds for the normal approximation of multiple integrals with respect to an isonormal 
Gaussian process. The bounds in [13] are also expressed in terms of contractions of the 
underlying kernel. 

Example 4.4 (Double integrals). We consider a double integral of the type 12(f), where 
/ E L 2 (/i 2 ) satisfies (|4.9p (according to Remark [4. 1( 1). this implies that (|4.8|) is satisfied). We 
suppose that the following three conditions are satisfied: (a) E^/) 2 = 2||/|| 2 2 ( M 2) = 1, (b) 
/ / G L 2 (f/) and (c) / G £ 4 (/x 2 ). Since q = 2 here, one has f(z, ■) *° /(z, -)(o) = f(z, a) 2 
which is square-integrable, and hence Assumption (ii) and condition (14, 16ft in Theorem 14.21 are 
verified. Using relations ()4.19p - (j4.20|) . we can deduce the following bound on the Wasserstein 
distance between the law of 12(f) and the law of X ~ ^(0, 1): 

d w (I 2 (f),X) (4.21) 

< Vs\\f*\ f\\ L *w + {2 + Vs(i + Ve)} \\f4 fWvw + ^f z j 4 d^. 

To obtain (14.2ip . observe first that, since assumption (a) above is in order, then relations 
(l4~T2l - (l4~T3l in the statement of Theorem gj yield that 

d w (h(f),x) < + gUD + 2 3 / 2 x (jHHD. 
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To conclude observe that 

6351) = 2{2 1 /2||/*i f\\ LW + 11/ f\\ L2M } 

and 

M) = 11/ 4 f\\m»){u 1/2 + 1} + 2 1 / 2 ||/ni 4(/ , 2) , 

since \\f *\ /Hl 2 (^) = 11/ *i /IIl 2 ^ 3 )- A general statement, involving random variables of the 
type F = I\(g) + h(h) is given in Theorem 16.11 n 



Example 4.5 (Triple integrals). We consider a random variable of the type 13(f), where 
/ G L^(fi 3 ) verifies (14. 8 1) (for instance, according to Remark l4.1f 4). we may assume that / 
has support contained in some rectangle of finite ^-measure) . We shall also suppose that 
the following three conditions are satisfied: (a) E/3(/) 2 = 3!||/||^w 3) = 1, (b) For every 
r = 1, 3, and every I = 1, rA2, one has that /■*£/£ L 2 (/x 6_r ~'), and (c) / G L 4 (/z 3 ). One 
can check that all the assumptions in the statement of Theorem 14.21 (in the case q = 3) are 
satisfied. In view of ()4.19|) - ()4.2Q|) . we therefore deduce (exactly as in the previous example and 
after some tedious bookkeeping!) the following bound on the Wasserstein distance between 
the law of I 3 (f) and the law of X ~ jV(Q, 1): 

d w (h(f),X) < 3(4! 1 /2)|| / ^ / || L2(M4) + (6 + 2 yi8)(3! 1 /2)|| / 4 / || i2{/j 3 ) (4.22) 

+6V2H/ 4 fhw) + {6 + Vl8(4 + 2(4! 1 /2))}|| / 4 /|| L2(/x2) 
+{v / 18(2 1 /2 + 5! 1 /2 )+6} || /+ 2 /||l2( ^ ) 

+{3 3 / 2 (18 1 / 2 )}^ 3 /^3. 

□ 



Proof of Theorem 14.21 First observe that, according to Theorem 13.11 we have that 

+ - / E\D z I q (f)\ 3 fi(dz). 



E 



--W D m\\h M 



d w (I q (f),X) < 
The rest of the proof is divided in four steps: 



1 



(SI) Proof of the fact that WE 
appearing at the line fj4.12H . 



(i-|ll^(/)lli 2(M) ) : 



is less or equal to the quantity 



(S2) Proof of the fact that - f z E\D z I q (f)\ s u.(dz) is less or equal to the quantity appearing 
at the line (l4~T3j) . 



(53) Proof of the estimate displayed in formulae (|4.14|) - (|4.15|) . 

(54) Proof of the estimate in formulae (|4.17p - (|4.18|) . under the assumption (|4.16|) . 
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Step (SI). Start by defining Af as the collection of those z G Z such that Assumption (ii) 
in the statement is violated, that is, 

A f = [zeZ:3r = l, q - 1, I = 0, r - 1 : f(z, •) *< f(z, •) g L 2 (/z 2M ~ r ~')} ■ 

By assumption, one has that fJ,(Af) = 0. By writing 

A\ \ = A f x • • • x A f G Z*, 

V „ ' 

q times 

one also deduces that u. q (A Q f ) = 0. As a consequence, I g (f) = I g (fl^), a.s.-P. It follows 
that, by replacing / with fl^i (and without loss of generality), we can assume for the rest of 
the proof that the following stronger assumption is verified: 

(if) For every z G Z, every r = 1, q — 1 and every I = 0, r — 1, one has that f(z, •) * l r 

fiz^TL^-V-r- 1 ). 



Now use (|2.17l) to write D z I q (f) = ql q ^i(f(z, •)), z G Z. By virtue of the multiplication 
formula (|2.14p . in particular (14, 2 1) . and by adopting the notation (14, ip . we infer that, for every 

ze z 



2(g-l) 



{D z I q (f)} 2 = qql [ f 2 (z, -)d^ + q 2 ]T I p (Gf l f(z, ■)). 
Jzq - 1 P =l 

Since (|4.8|) is in order, one has that, for every p = 1, 2(g — 1), 

E j z \I p {Gl- l f{z,-))\u.(dz) ^ Up\ J^G q p - 1 f(z,-) 2 duPu.(dz) < oo, 



(4.23) 



(4.24) 



where we have used the Cauchy-Schwarz inequality, combined with the isometric properties 
of multiple integrals, as well as the fact that, for every z G Z, 

Gf^fiz,- 1 



Zp 



V< [ GfV&OV 
Jzp 



Relations (|4.23p - (l4.24l) yield that one can write 



\DI q (f)\\h (p) 



- [ {D z I q (f)} 2 ^(dz) - 1 
1 J z 

2(3-1) . 

1 + 9 E / Ip{Gr X f{z,-))^dz). 
p=l Jz 



(4.25) 



Since assumption (14.8(1 is in order, one has that, for every p = 1, 2(g — 1), 



E 



i p (G*-7(s,0M<fc; 



(4.26) 



< y E |/ P (G«" 1 /(^, -))I p {Gl- l f{z', -))| n(dz)u.(dz') 
^pl UJ J^i^ttzrWdtjP »(dz)\ <oo, 



19 



and one can easily verify that, for 1 ^ p ^ I ^ 2(q — 1), the random variables 

I p (Gf l f{z, -))fi(dz) and / T^G?" 1 f(z, -))p(dz) 
z Jz 

are orthogonal in L 2 (P). It follows that 

2 



E 



1 



1 



DI q (f)\\h M 
so that the estimate 



2(5-1) 

mf\\h { , q) -i) 2 +Q 2 E E 

p=i 



i^g^/Ov))/^) 



E 



i-Jn^(/)iiW))' 



< (mi 



is proved, once we show that 

f /-(Gg" 1 /^,-))/*^) 



{d*f} 2 dv?. 



(4.27) 



The proof of (|4.27j) can be achieved by using the following relations: 



E 



Z 



E[I p (G£-7(^ ■)I P (Gf 1 f(z', -)Mdz)^dz') 



(4.28) 



z Jz 



ZP 



GrV(^-)Gr7(^,-)^ p 



g-l 



Gl~ X f{z,-)n(dz) 



dn p ^ p\ 



fj,(dz)p,(dz') 

Gf x f(z r )p,{dz) 



dfjf 



Zp 



{G*/}V r 



Note that the use of the Fubini theorem in the equality (|4.28|) is justified by the chain of 
inequalities (|4.26|) . which is in turn a consequence of assumption (14.81) . 
Step (S2). First recall that 



E 



Now use the Cauchy-Schwarz inequality, in order to write 



q 



\D z I q (f)\ 3 fi(dz) 



^ -WE 



\DI q (f)\\h 



(p) 



I j E[(D 9 I q {f))*]fi(dz) (4.29) 



(9-1)111/11^^) x J Jn{DJ q {f)fUdz) 



(4.30) 



By using (|4.4|) and (I4.23p . one deduces immediately that 

2(g-i) 

{D z I q (f)f = q 2 £ JplGTVM. 

p=0 
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As a consequence, 



2(g-i) 



1 jnDMfYUdz) < q 2 Jpljfj^ 1 



\ P=0 
2(9-1) 

< ? 2 E ^ 1/2 < 

p=0 



(4.31) 



Z VJZv 



yielding the desired inequality. 

Step (S3). By using several times the inequality ya + b ^ 
in order to prove (I4,14j) - (|4.15p . it is sufficient to show that 

2(9-1) j—. : 

1} <x-t-.w,-iyw($-_l)£-_\ 



a, b ^ 0) one sees that, 



q iA(g-l) 



t=l s=l 



1/2 



To see this, use (|4.7p and the fact that (by (|2.6|) ) ||/*f /Hl 2 ^ 2 ?-*-*) ^ ||/*| /llz, 2 ^ 2 ?-*-")) to 
obtain that 



Lxp 



1 2 
/} ^ 



j II/*? /llL 2 (/i 2 9-t-s), 



t-l/Vs-1 



2(9-1) 

E ^\ 

P =i 

2(9-1) g tA(g-l) 

< E p'- 1/2 E E i {2 9-^ } (*-i)! 1/2 

p=l t=l s =l 

and then exploit the relation 

2(g-l) g tA(q-l) q tA(g-l) 

E ^ /2 E E i{29-t- s = P} = E E w<29-i}(2.-*- s )! 1/2 - 

p=l t=l s =l t=l s=l 

Step (S4). By using (14. 3h and some standard estimates, we deduce that 



2(9-1) 

E 

P =o V ^ ^ 

2(9-1) 9-1 r 

< E P !l/2 EE 1 {2(9-l)-r- i = P }^ 
p=0 r=0 Z=0 



,1/2 



1/2 



(V)C) * 



(4.32) 



[/(z.OH/C^OF^K^)- ( 4 - 33 ) 



Z JZP 



We claim that, if (|4.16p is satisfied, then, for every r = 0, ...,q — 1 and Z = 0, ...,r 

[f( z ,.)4f(z,-)] 2 d^^dz) = [ [f* 9 q z]~ r fW +r+1 . 

Jz'+ r + 1 




(4.34) 



'Z JZP 

In the two 'easy' cases 



21 



(a) r = q — 1 and I = 1, q — 1, 



(b) r = 1, g — 1 and Z = 0, 

relation (j4.34j) can be deduced by a standard use of the Fubini theorem (assumption (|4.16p is 
not needed here). Now fix p = 1, 2q — 2, as well as r = 1, ...q — 2 and I = 0, ...,r in such a 
way that 2(q — 1) — r — I = p. For every fixed z G Z, write \ f(z, -)| * l r \ f(z, -)| to indicate the 
contraction of indices (r, I) obtained from the positive kernel \ f(z, -)|. Note that, for z fixed, 
such a contraction is a function on Z p , and also, in general, \ f(z, -)\-k l r \f(z, -)| ^ f(z, -)~k l r f(z, •) 
and -)| *r 1/(^1 ')l 7^ f( z > ')*r f{ z i ')• By a standard use of the Cauchy-Schwarz inequality 
and of the Fubini theorem, one sees that 




|/0V)|4I/0V)I d^n(dz) ^ ||/*r 1 - r+, /||^ r - I -i ) <oo, (4.35) 



.q-l-r+l f\\2 

where the last relation is a consequence of assumption (14. 16[) as well as of the fact that, by 
construction, 1^1 + r — I ^ q — 1. Relation (|4.35j) implies that one can apply the Fubini 
Theorem to the quantity 

/ / [f(z r )* l r f(z r )] 2 d^(dz) 

JZ JZP 

(by first writing the contractions f(z, •) * r f(z, •) in an explicit form), so to obtain the desired 
equality (14. 34p . By plugging (14, 34p into (|4.32l) and by applying the change of variables a = 
q—l — r and b = q — l, one deduces (|4.17p - (|4.18p . This concludes the proof of Theorem 14.21 □ 



5 Central limit theorems 

We consider here CLTs. The results of this section generalize the main findings of [23]. The 
following result uses Theorem 14.21 in order to establish a general CLT for multiple integrals of 
arbitrary order. 

Theorem 5.1 (CLTs on a fixed chaos) Let X ~ ^(0,1). Suppose that (J,(Z) = oo, fix 
q 2, and let Fk = I q (fk), k ^ 1, be a sequence of multiple stochastic Wiener-Ito integrals 
of order q. Suppose that, as k — > oo, the normalization condition ¥,(F 2 ) = ^H/fcll^^q) - ► 1 
takes place. Assume moreover that the following three conditions hold: 

(I) For every k ^ 1, the kernel fk verifies (4-^ f or every p = 1, ...,2(q — 1). 

(II) For every r = 1, ...,q, and every I = 1, ...,rA(q — 1), one has that fk* l r fk£ L 2 (fi 2q ~ r ~ l ) 
and also Wfk^i fk\\L 2 {^ 2 i- r — 1 ) ~^ ( as k ^ oo). 

(Ill) For every k ^ 1, one has that J Zq f£ dfi q < oo and, as k — > oo, f Zq f£ dpfl — ► 0. 



Then, Fk ^+ X , as k — > oo, and formulae (f^. 14 )- \4.20^ provide explicit bounds in the Wasser- 
stein distance dw(Iq(fk), X). 

Proof. First note that the fact that * r fk G L 2 (// 2?_ ) for every r = 1, ...,q and every 
I = 1, ...,r A (q — 1) (due to Assumption (II)) imply that Assumption (ii) in the statement 
of Theorem 14.21 holds for every k ^ 1 (with fk replacing /), and also (by using (|4.19l) ) that 
condition (|4.16p is satisfied by each kernel fk- Now observe that, if Assumptions (I)- (III) are in 
order and if E(F fe 2 ) = q\\\fk\\ 2 L 2 i/iq) -> 1, then relations (f4Tl4) - (lO0l) imply that d w (F k ,X) -> 0. 
Since convergence in the Wasserstein distance implies convergence in law, the conclusion is 
immediately deduced. 
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□ 



Example 5.2 Consider a sequence of double integrals of the type h(fk), k ^ 1 where fk £ 
Lg(n ) satisfies (|4.9p (according to Remark 14.11 (1). this implies that (|4.8|l is satisfied). We 
suppose that the following three conditions are satisfied: (a) E/2(/fc) 2 = 2||/fc||j : 2( A( 2) = 1, (b) 
/fc*2 fk ^ -^ 2 (m 1 ) an( l ( c ) fk ^ L 4 (/j?). Then, according to Theorem 15.11 a sufficient condition 
in order to have that (as k — > oo) 

/ 2 (/ fc )-^(0,l) (5.1) 
is that 

ll/fcll^) - 0, (5.2) 

ll/fc 4 /fclU 2 0) -> 0, and ||/ fc *] /fe||i2 (At 2) -> 0. (5.3) 

This last fact coincides with the content of Part 1 of Theorem 2 in [23], where one can find 
an alternate proof based on a decoupling technique, known as the "principle of conditioning" 
(see e.g. Xue (34]). Note that an explicit upper bound for the Wasserstein distance can be 
deduced from relation A4.21 j) . □ 



Example 5.3 Consider a sequence of random variables of the type hifk), k ^ 1, with unitary 
variance and verifying Assumption (I) in the statement of Theorem 15.11 Then, according to 
the conclusion of Theorem 15.11 a sufficient condition in order to have that, as k — > oo, 

h(f k ) ^ x ~ jr(0A), 

is that the following six quantities converge to zero: 

ll/fc *1 fk\\z*(ji*), ll/fc *2 /fc|U 2 Gu 3 )> ll/fc *2 fk\\L 2 (fi 2 ), 
ll/fc *3 /fcllL 2 (At 2 )' ll/fc *3 /fclli 2 ^ 1 )' and ll/fclli 4 ^ 3 )- 

Moreover, an explicit upper bound in the Wasserstein distance is given by the estimate (|4.22|) . 

□ 



The following result, proved in [23j Theorem 2], represents a counterpart to the CLTs for 
double integrals discussed in Example 15.21 

Proposition 5.4 (See [23]) Consider a sequence = hifk), k ^ 1, of double integrals 
verifying assumptions (a), (b) and (c) of Example \5. 6 A Suppose moreover that 15. ty) takes 
place. Then, 

1. if E L 4 (P) for every k, a sufficient condition to have 115. 3\) is that 

E(F 4 )^3; (5.4) 

2. if the sequence {F^ : k ^ l} is uniformly integrable, then conditions h5. 115. 3\) and 
J5.^p are equivalent. 

For the time being, it seems quite hard to prove a result analogous to Proposition 15.41 for 
a sequence of multiple integrals of order q ^ 3. 
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6 Sum of a single and a double integral 



As we will see in the forthcoming SectionEJ when dealing with quadratic functionals of stochas- 
tic processes built from completely random measures, one needs explicit bounds for random 
variables of the type F = I\(g) + 12(h), that is, random variables that are the sum of a single 
and a double integral. The following result, that can be seen as a generalization of Part B of 
Theorem 3 in [23], provides explicit bounds for random variables of this type. 

Theorem 6.1 Let F = h(g) + h{h) be such that 
(I) The function g belongs to L 2 (fi) n L 3 (/i); 

(II) The kernel h G L 2 (fi 2 ) is such that: (a) h*\h G L 2 (fi 1 ), (b) relation |^.g[ ) is verified, 
with h replacing f, and (c) h G L 4 (/x 2 ). 

Then, one has the following upper bound on the Wasserstein distance between the law of F 
and the law of X ~ ^(0, 1): 



dw(F,X) 



(6.1) 



1 



+ 2\\h*l h\\ L 2 



(/■<) 



(6.2) 



MI£ 2(M) -2| 

+V8\\h 4 h\\ L 2 { ^ + % 4 h\\ L 2 (fl) + 32\\gf L3M 
+4117*11^) x {\\h\\ 2 LW + 2 1 ^\\h4h\\ L2(p) }. 
The following inequality also holds: 

\\g*\ ^IIl 2 ( M ) < llslU 2 ao x Hl^y 

Proof. Thanks to Theorem 13.11 we know that d\y(F, X) is less or equal to the RHS of (13.21) . 
We also know that 

D Z F = g(z) + 2I x {h(z, •)), and — D z L~ l F = g(z) + h(h(z, •)). 

By using the multiplication formula (|2.14p (in the case p = q = 1) as well as a Fubini argument, 
one easily deduces that 

j D Z F x (-D z L~ l F)ii{dz) = \\g\\l^ ) +2\\h\\l HlJ?) +2I 1 (h-k\ h) + 2I 2 {h*\ h)+$h{g*\ h). 
J z 

This last relation yields 



E[(l-(DF,-DL-^F) L 2 M y 



(6.3) 



\9\\ 2 L 2 itl) -2\\h» 2 



L2( M 2) 



+ 2\\h*l h\\ L 2 ip) + y/8\\h*\ /i|| L 2 (At2) + 3\\g*{ h\\ L 2 M . 
To conclude the proof, one shall use the following relations, holding for every real a, b: 

(a + 26) 2 |a + b\ < (|a| + 2|6|) 2 (|a| + |6|) < (\a\ + 2\b\f ^ 4|a| 3 + 32|6| 3 . (6.4) 
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By applying (|6.4p in the case a = g(z) and b = I±(h(z, •)), one deduces that 

E [\D Z F\ 2 \D Z L- 1 F\] fi(dz) (6.5) 



z 



^4E / \h(h(z,-))\^(dz) + 32 / \g(z)\ 6 fi(dz). 
Jz Jz 

By using the Cauchy-Schwarz inequality, one infers that 

4E / \h(h(z,-))\ 3 fi(dz) 
Jz 



< WE / \h(h(z,Wn(dz) x 4|H| W) , (6.6) 



z 

and (16. ip is deduced from the equality 

z 



Formula (|6.2D is once again an elementary consequence of the Cauchy-Schwarz inequality. □ 



Remark 6.2 Consider a sequence of vectors {F k ,H k ), k ^ 1 such that: (i) Fj. = h{fk), 
k ^ 1, is a sequence of double integrals verifying assumptions (a), (b) and (c) in Example 15.21 
and (ii) H k = Ii(hk), k ^ 1, is a sequence of single integrals with unitary variance. Suppose 
moreover that the asymptotic relations in (|3.7|) . (|5.2p and (|5.3|) take place. Then, the estimates 
(I6.ip - (|6.2p yield that, for every (a, (3) ^ (0,0), the Wasserstein distance between the law of 

J— (aF k + (3H k ) 
\Ja 2 + f3 2 

and the law of X ~ ^(0, 1), converges to zero as k — > oo, thus implying that (F^,i/fc) 
converges in law to a vector (X, X') of i.i.d. standard Gaussian random variables. Roughly 
speaking, this last fact implies that, when assessing the asymptotic joint Gaussianity of a 
vector such as (Fk,Hk), one can study separately the one-dimensional sequences {F k } and 
{H k }. This phenomenon coincides with the content of Part B of Theorem 3 in [23J. See |15| . 
[T7] and [21] for similar results involving vectors of multiple integrals (of arbitrary order) with 
respect to Gaussian random measures. 

7 Applications to non-linear functionals of Ornstein-Uhlenbeck 
Levy processes 

As an illustration, in this section we focus on CLTs related to Ornstein-Uhlenbeck Levy pro- 
cesses, that is, processes obtained by integrating an exponential kernel of the type 



2Axe- A ^l Kt} , 

with respect to an independently scattered random measure. Ornstein-Uhlenbeck Levy pro- 
cesses have been recently applied to a variety of frameworks, such as finance (where they are 
used to model stochastic volatility - see e.g. p]) or non-parametric Bayesian survival analysis 
(where they represent random hazard rates - see e.g. [4j [12j [20]). In particular, in the refer- 
ences [1] and [2D] it is shown that one can use some of the CLTs of this section in the context 
of Bayesian prior specification. 
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7.1 Quadratic functionals of Ornstein-Uhlenbeck processes 

We consider the stationary Ornstein-Uhlenbeck Levy process given by 



Y t A = V2X uexp(-\(t-x))N(du,dx) , t ^ 0, (7.1) 

J \ Jm. 

where N is a centered Poisson measure over Ixl, with control measure given by v (du) dx, 
where v {■) is positive, non-atomic and <7-finite. We assume also that f u J v(du) < oo for 
j = 2,3,4,6, and J u 2 v{du) = 1. In particular, these assumptions yield that 



E 



yf 



2A / / u 2 e~ 2X[t ~ x ^v (du) dx = 1. 



The following result has been proved in [23j Theorem 5]. 
Theorem 7.1 (See |23j) For every A > 0, as T — > oo, 



Q(T,X):=VT^£ (Y*) 2 dt - 1 } *Z> ^\+cl x X, (7.2) 

where c 2 , := J u 4 v(du) and X ~ jV (0, 1) is a centered standard Gaussian random variable. 

By using Theorem 16.11 one can obtain the following Berry-Esseen estimate on the CLT 
appearing in (|7.2)) (compare also with Example 13 .6p . 



Theorem 7.2 Let Q(T,X), T > 0, be defined as in VTE) . and set 



Q(T,X) :=Q(T,A)/^i + C 2. 
Then, there exists a constant < 7(A) < 00, independent ofT and such that 

d w (Q(T,X),X)^^. (7.3) 



(7.4) 



Proof. Start by introducing the notation 



H X>T (u, x; u,' x>) = (u x „') '^!l X0 {e^+*') (l - e~ 2T ) l (xVx%0) + 



Tj\ + c 2 



, Mx+x') f-2\(xVx') C -2XT\ -, \ 

+ e ' [e - e J 1( X \/ X < >0 )J , 

tt* (n, rA — „,2 ~*-(-°o,T] ( X ) f „2\x (-i C ~2T\ -, 1 „2\x f —2Xx -2\T\ \ 

H XT {u,x) - u je (1-e J 1(^0) + e -e )l( x>0 )j. 



TJi + c 2 



As a consequence of the multiplication formula (I2.14p and of a standard Fubini argument, one 
has (see [23j Proof of Theorem 5]) 

Q{T, A) = h (VtHIt) + h (VTh^t) , 



a combination of a single and of a double integral. To apply Theorem 16. 11 we use the following 
asymptotic relations (that one can verify by resorting to the explicit definitions of H\^p and 
H\ t gi ven m |7.4p ). holding for T — > 00: 
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(a) 



(b) 



(d) 



(f) 



Vth$ 



L 2 (dudx) 



Vth 



Vth 



L 2 ({dvdx) 2 ) 



L 3 (dudx) y/T 



Vth x 



T 



^{(dvdx)) 2 y/T 

[Vth x>t ) 4 (VfH XyT ) 



(Vth x , t ) *} (Vth x>t ) 



L 2 {dvdx) y/T 



L 2 {(dvdx) 2 ) 



THt T ) *} (VTH XsT ] 



L 2 {dvdx) y/T 

The conclusion is obtained by using the estimates (|6.ip - (l6.2l) and applying Theorem 16.11 ^ 



7.2 Berry-Esseen bounds for arbitrary tensor powers of Ornstein-Uhlenbeck 
kernels 

Let ./V be a centered Poisson measure over KxK, with control measure given by fi(du, dx) = 
v (du) dx, where v (•) is positive, non-atomic and cr-finite. We assume that f u 2 v (du) = 1 and 
/ u 4 v (du) < oo. Fix A > 0, and, for every t ^ 0, define the Ornstein-Uhlenbeck kernel 

ft(u,x) = u x V2Xexp{-X(t - x)}l{ x ^t}, (it, x) e R x R. (7.5) 

For every fixed q ^ 2, we define the gth tensor power of /t, denoted by ff q , as the symmetric 
kernel on (R x R) 9 given by 



ff q (u 1 ,x 1 ;...;ug,x q ) = Y[ft(uj,Xj). (7.6) 

We sometimes set y = (u,x). Note that, for every t ^ 0, one has that J ft(y)fi(dy) = 
J ft(u,x)v(du)dx = 1, and therefore ff q G L 2 (/j, q ); it follows that the multiple integral 

Zi(g) := / ff (/f ? ) 

is well defined for every t ^ 0. 
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Remark 7.3 Fix q ^ 2, and suppose that J \u\idv < oo, Vj = 1, 2q. Then, one can prove 
that Z t (l) q is square-integrable and also that the random variable Z t (q) coincides with the 
projection of Z t (l) q on the qth Wiener chaos associated with N. This fact can be easily 
checked when q = 2: indeed (using Proposition 12.51 in the case f = g = ft) one has that 

Z t {\? = h(ff) + J 2 (/f 2 ) = h{ft) + Z t (2), 

thus implying the desired relation. The general case can be proved by induction on q. 

The main result of this section is the following application of Theorem 14.21 and Theorem 

Theorem 7.4 Fix A > and q ^ 2, and define the positive constant c = c(q,X) := 2(q — 1)!/A. 
Then, one has that, as T — > oo, 

rp 

1 / / \ t, law 



M T (q) := — I Z t (q)dt ^ X ~ JT$, 1), (7.7) 



n 



and there exists a finite constant p = p(\, q,v) > such that, for every T > 0, 

d w (M T (q), X) < P (7.8) 



Proof. The crucial fact is that, for each T, the random variable M^(g) has the form of a 
multiple integral, that is, Mx(q) = I q (FT), where Ft S L 2 s {p q ) is given by 



F T (ui,xr,...;u q ,x q ) = —= \ f t q (u u xr, u q , x q )dt, 



T 





where ff q has been defined in (|7.6p . By using the fact that the support of Ft is contained 
in the set (Ix (— oo,T]) q as well as the assumptions on the second and fourth moments of 
u, one easily deduces that the technical condition (|4.8p (with Ft replacing /) is satisfied for 
every T ^ 0. According to Theorem 14.21 and Theorem 15.11 both claims (|7.7p and (|7.8|) are 
proved, once we show that, as T — > oo, one has that 

|1 -E{M T (q)) 2 \ ~ 1/T, (7.9) 

and also that 

\\FT\\h^) = 0(l/n and (7.10) 

\\F T * l r F T \\i? 0j ? 9 -r-i ) =O(l/T), Vr = l, VI = l,..,rA( ? -l) (7.11) 

(relation (|7.9|) improves the bounds in Theorem 14. 2p . In order to prove (|7.9p - (l7.1ip . for every 
ti,t2 ^ we introduce the notation 

(/t 1 ,/t 9 >M= / ftAy)My)Kdy) = e- x ^ +t ^e 2X ^ At *\ (7.12) 

JMxE 

(recall that f u 2 dv = 1) and also, for £1,^2, £3^4 ^ 0, 

(ftj t2 , f t Ju)» = I ftMftMftMftM Kdy) (7.13) 



u 4 u(du) 



x ^ e -A(ti+t 2 +i3+t4) e 4A(tiAi 2 At3Ai4) ^7 
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To prove (|7.9|) . one uses the relation (|7.12p to get 



E [M T {qf 



I /-T /-T 



cT 



^0 



(ft 1 ,ft 2 )ldt 1 dt 2 = l 



TqX 



(1-e 



In the remaining of the proof, we will write k in order to indicate a strictly positive finite 
constant independent of T, that may change from line to line. To prove (|7.10j) . one uses the 
fact that 



t 

[ F T d^ = ^[ [ [ (ftJt^ftMldhd^dtsdU^^, 
J(txl)i c 1 Jo Jo Jo Jo 1 

where the last relation is obtained by resorting to the explicit representation (|7.14p . and then 
by evaluating the restriction of the quadruple integral to each simplex of the type {^(i) > 
^tt(2) > *tt(3) > *7r(4)}> where ir is a permutation of the set {1,2,3,4}. We shall now verify the 
class of asymptotic relations (|7.11|) for r = q and I = 1, q. With y = (u, x), one has 



T r T 



q-l 

Y[fti {vi)ft 2 {yi 

,i=l 



(ftjt 2 )^ dhdt 2 , 



and hence 



\\F T * l q F T f m 



K 
K 



q * J \\L*(p,<l- 1 ) 
T r T rT r T 



JO JO JO 



(ftjt 2 , ftJu)V(ftJt 2 )l(ft 3 ftX dhdhdhdU 



where the last relation is verified by first using (|7.12p - (l7,14j) . and then by assessing the re- 
striction of the quadruple integral to each one of the 4! = 24 simplexes of the type {^(i) > 
^tt(2) > *7r(3) > *7r(4)}- To deal with (17. lip in the case r = 1, q — 1 and / = 1, r, one uses 
the fact that 



Ft * r Ft(vi, ...,y r -i,wi, ...,w q - r ,zi, ...,z t 



q—r ) 



1 

cT 



T rT 
JO 



r-l 



\ft 1 {yi)ft 2 {Vi) 



i=l 



X (fhiwi) • ■ • f tl {w q - r )ft 2 {zi) • • • ft 2 (z q - r )) {ftjt 2 )i dtidt 2l 



and therefore 



\\F T * r F T \\ 2 L2 ^ 2q - r -l) 



K 

T*J 

K 

T 



T rT rT rT 



JO JO 



{ftj t2 , /* 3 / t4 )r i (/*i/t3)r r (/* 2 /*4)r r (Ai/t 2 )j t (/t3/ i 4)M dhdt 2 dt 3 du 



where the last relation is once again obtained by separately evaluating each restriction of the 
quadruple integral over a given simplex. This concludes the proof. ^ 



29 



Acknowledgments. Part of this paper has been written while G. Peccati was visiting the 
Department of Mathematics of the "Universitat Autonoma de Barcelona", in April 2008. This 
author heartily thanks J.L. Sole and F. Utzet for their hospitality and support. J.L. Sole and 
F. Utzet were supported by the Grant BFM2006-06427, Ministerio de Education y Ciencia 
and FEDER. M.S. Taqqu was partially supported by the NSF Grants DMS-0505747 and 
DMS-0706786 at Boston University. 

References 

[1] Bandorff-Nielsen, O.E. and Shepard, N. (2001). Non Gaussian Ornstein-Uhlenbeck-based 
models and some of their uses in financial economics. Journal of the Royal Statistical 
Society B 63, 167-241. 

[2] Chatterjee, S. (2007). A new method of normal approximation. To appear in: The Annals 
of Probability. 

[3] Chen, L. and Shao, Q.-M. (2005). Stein's method for normal approximation. In: An 
introduction to Stein's method, 1-59. Lect. Notes Ser. Inst. Math. Sci. Natl. Univ. Singap. 
4, Singapore Univ. Press, Singapore, 2005. 

[4] De Blasi, P., Peccati, G. and Priinster, I. (2008). Asymptotics for posterior hazards. To 
appear in: The Annals of Statistics. 

[5] Decreusefond, L. and Savy, N. (2007). Rubinstein distance on configuration space. 
Preprint. 

[6] Dudley, R.M. (2002). Real Analysis and Probability. Cambridge University Press. Cam- 
bridge. 

[7] ltd, K. (1956). Spectral type of the shift transformation of differential processes with 
stationary increments. Trans, of the American Mathematical Society 81, 252-263 

[8] Kabanov, Y. (1975). On extended stochastic integrals. Theory of Probability and its Ap- 
plications 20, 710-722. 

[9] Kyprianou, A.E. (2006). Introductory Lectures on Fluctuations of Levy Processes with 
Applications. Springer- Verlag. Berlin Heidelberg New York. 

[10] Lee, Y.-J. and Shih, H.-H. (2004). The product formula of multiple Levy-Ito integrals. 
Bull. Inst. Math. Acad. Sinica 32, 71-95. 

[11] J. Leon, J.L. Sole, F. Utzet and J. Vives (2002). On Levy processes, Malliavin calculus 
and market models with jumps. Finance and Stochastics 6, 197-225. 

[12] Nieto-Barajas, L.E. and Walker, S.G. (2005). A semi-parametric Bayesian analysis of 
survival data based on Levy-driven processes. Lifetime Data Analysis 11, 529-543. 

[13] Nourdin, I. and Peccati, G. (2008a). Stein's method on Wiener chaos. To appear in: 
Probability Theory and Related Fields. 



30 



[14] Nourdin, I. and Peccati, G. (2008b). Stein's method and exact Berry-Esseen bounds for 
functionals of Gaussian fields. Preprint. 

[15] Nourdin, I., Peccati, G. and Reveillac, A. (2008). Multivariate normal approximation 
using Stein's method and Malliavin calculus. Preprint. 

[16] Nualart, D. (2006). The Malliavin calculus and related topics. 2 nd Edition. Springer- 
Verlag. Berlin Heidelberg New York. 

[17] Nualart, D. and Ortiz-Latorre, S. (2007). Central limit theorems for multiple stochastic 
integrals and Malliavin calculus. Stochastic Processes and their Applications 118(4), 614- 
628. 

[18] Nualart, D. and Peccati, G. (2005). Central limit theorems for sequences of multiple 
stochastic integrals. The Annals of Probability 33 (1), 177-193. 

[19] Nualart, D. and Vives, J. (1990). Anticipative calculus for the Poisson process based on 
the Fock space . Seminaire de Probabilites XXIV, LNM 1426, 154-165, Springer, Berlin. 

[20] Peccati, G. and Priinster, I. (2008). Linear and quadratic functionals of random hazard 
rates: an asymptotic analysis. To appear in: The Annals of Applied Probability. 

[21] Peccati, G. and Taqqu, M.S. (2007). Stable convergence of I? generalized stochastic 
integrals and the principle of conditioning. The Electronic Journal of Probability, 12, 
447-480, n. 15 (electronic). 

[22] Peccati, G. and Taqqu, M.S. (2007). Limit theorems for multiple integrals. Preprint. 

[23] Peccati, G. and Taqqu, M.S. (2008). Central limit theorems for double Poisson integrals. 
To appear in: Bernoulli 

[24] Peccati, G. and Tudor, C.A. (2004). Gaussian limits for vector- valued multiple stochastic 
integrals. In: Seminaire de Probabilites XXXVIII, 247-262, Springer Verlag 

[25] Picard, J. (1996). Formules de dualites sur l'espace de Poisson. Ann. Inst. H. Poincare 
B, 32(4), 509-548. 

[26] Reinert, G. (2005). Three general approaches to Stein's method. In: An introduction to 
Stein's method, 183-221. Lect. Notes Ser. Inst. Math. Sci. Natl. Univ. Singap. 4, Singapore 
Univ. Press, Singapore. 

[27] Sato, K.-I. (1999). Levy Processes and Infinitely Divisible Distributions. Cambridge Stud- 
ies in Advanced Mathematics 68. Cambridge University Press. 

[28] Sole, J.L., Utzet, F. and Vives, J. (2007). Canonical Levy processes and Malliavin calcu- 
lus. Stochastic Processes and their Applications 117, 165-187. 

[29] Stein, Ch. (1972). A bound for the error in the normal approximation to the distribution of 
a sum of dependent random variables. In: Proceedings of the Sixth Berkeley Symposium 
on Mathematical Statistics and Probability, Vol. II: Probability theory, 583-602. Univ. 
California Press, Berkeley, Calif.. 



31 



[30] Stein, Ch. (1986). Approximate computation of expectations. Institute of Mathematical 
Statistics Lecture Notes - Monograph Series, 7. Institute of Mathematical Statistics, 
Hay ward, CA. 

[31] Surgailis, D. (1984). On multiple Poisson stochastic integrals and associated Markov 
semigroups. Probabilty and Mathematical Statistics 3(2), 217-239. 

[32] Surgailis, D. (2000). CLTs for Polynomials of Linear Sequences: Diagram Formulae with 
Applications. In: Long Range Dependence, 111-128, Birkhauser. 

[33] Surgailis, D. (2000). Non-CLT's: U-Statistics, Multinomial Formula and Approximations 
of Multiple Wiener-Ito integrals. In: Long Range Dependence, 129-142, Birkhauser. 

[34] Xue, X.-H. (1991). On the principle of conditioning and convergence to mixtures of dis- 
tributions for sums of dependent random variables. Stochastic Processes and their Appli- 
cations 37(2), 175-186. 



32 



